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Abstract 

A method is presented to calculate the spin relaxation times Ti , T2 due to a non- 
uniform magnetic field, and the linear-in-electric-field precession frequency shift 
6uje when an electric field is present, in the diffusion approximation for spins 
confined to a rectangular cell, ft is found that the rectangular cell geometry 
admits of a general result for Ti, T2, and 5uje in terms of the spatial cosine- 
transform components of the magnetic field. 
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1. Introduction 

An experiment to measure the neutron electric dipole moment (nEDM), 
to be installed at the FnPB beamline at Oak Ridge National Laboratory, will 
utilize a helium-3 comagnetometer in the central, superfluid-helium-filled mea- 
surement cell ^ [2] • The helium-3 polarization must remain high over the entire 
measurement period, ~1000 seconds, as the helium atoms precess in the hold- 
ing field and diffuse within a rectangular cell. Also, as there is a strong electric 
field E applied across the cell, a subtle effect, in which the interplay of the 
motional v x E field with gradients in the static magnetic field cause the pre- 
cession frequency to shift linearly with E 3J , must be well-understood or shown 
to be negligible. Design optimization of the experimental apparatus includes 
calculating the helium-3 spin relaxation times Ti , T2 , and linear-in-electric-field 
frequency shift 6uje due to given magnetic field non- uniformities. 

In this article, a method is shown to calculate these quantities in the dif- 
fusion approximation in a rectangular cell and for an arbitrary magnetic field. 
The starting point for the relaxation times is the Redfield theory of spin re- 
laxation 4 . In second order perturbation theory these can be written, for a 
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holding field in the z direction, as [S] 



7^ = 7^ (^2;i^('^o) + fcaa(wo)) , (1) 
;^ = ;^+7'fc..(0), (2) 

J 2 ^Jl 

where the spectral density is given in terms of magnetic field perturbations 

1 

kij{uj) ^ - / {hi{t)hj{t + t)) cos uirdr (3) 

^ J -oo 

Here, the total field in each direction i is Hi{t) — {Hi{t)) + hi{t) , such that aver- 
age perturbation {hi{t)) — 0, and wq = -f{Hz{t)) is the average spin precession 
frequency. 

McGregor [6j calculated the ensemble average correlation of the field pertur- 
bations seen by a diffusing particle in the case of a time-independent, uniform 
gradient of Hz in the x-direction, 

{hzit)Kit + T))^(^^^ {X{t)x{t + T)), (4) 

resulting in an analytic expression for T2 in a rectangular cell, 



T2 2Ti 120£i 



dx 



We relax the requirement of uniform gradient and find, in the case of a rect- 
angular prism cell, that Ti and T2 can be written in terms of the components 
of the 3D cosine transform of hq{r) over the cell volume. The same technique is 
applied to dressed spins [7], with uniform holding field and non- uniform dress- 
ing field, by mapping non-uniformities in the dressing field to equivalent non- 
uniformities in the holding field. In Section [3] a variation of the technique is 
used for the linear-in-electric-field frequency shift. Finally, as an example, in 
Section |4] the relaxation times and linear-in-electric-field frequency shift due to 
a superconducting rod near the cell are computed. 



2. Correlation functions in the diffusion limit 

The correlation function of hi can be expressed as integrals over the cell 
volume weighted by the probability density p{fo,t) that the particle is at fg at 
the initial time Iq, and the joint probability density p(f*, i|roj^o) that a particle 
at To at time to will be at f at time t. Thus, [6 

{hi{to)hi{tn + t)) = / dfnhi{ro)p{fo,to) / drh.i{f)p{r,tn + T\ro,to). (6) 
Jv Jv 
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The particle density will be taken as uniform in the cell, p{r,t) = 1/V. The 
joint probability is the solution to the diffusion equation, 

^p(f,i|ro,to) - DVMr,t\fo,to), (7) 

subject to reflecting boundary conditions at the walls, 

Vp(fs,t|ro,to)-n-0, (8) 

where rs & S and h is normal to the wall. In a rotated coordinate system 
(a;', y' , z') aligned with the cell walls, the solution for a x Ly x Lz box with 
walls at x' = y' = zLLy/2 and z' = zLLz/2 is 

p(r', t\r'o,to) = p{x', t\xQ, to; L.x)p{y', t\yQ, to; Ly)p{z' , to; L^), (9) 

with the ID solution dependent on the time difference t = t — to, [1] 

p{x',t\x'o,to;L,) = J^ + J^ e " " sm — sm 

+ — 2^ e " ^^z-^- cos — — cos— (10) 

^ n=2,4,... ^ 

It will be convenient to recognize the following: 

p{x'-L^/2, t\x'o~L,/2, to;L^) = ^ + ^ ^ e"" " ^"/^- cos cos — ^ 

n— 1,2,3,... 

(11) 

Putting Eq. [ojand p(rQ,to) — l/V into Eq. [6j changing the limits of integra- 
tion to < (/,■ < Li for each dimension and using Eq. II we have 



i(T) = {h,{to)h,{to + t)) 

OO 

= E 



OO 

e 

ria; ,ny — 

X dx^ dz^, - L.j,/2, y'o - Lj,/2, z'q - L J2) 

^ ^ ^ f^2;7^a;;, riyiry'^ n^nz'^ 

Cn^ Cny Cn, COS — — ^ COS ^ COS ^ 

-t^x -t^j/ -t^z 

X:^ / dx' dy' dz' hi{x' - L^/2,y' - Ly/2,z' - LJ2) 



V 



v 



r' n n '^^^^ ^v^V "^'^^ /'io\ 

C'nj, Cn, COS — COS COS — (12) 

J^y J-iz 



in which the factor C„ has been introduced. 



^" " I V2 otherwise. ^^^^ 
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We identify the 3D cosine transform of hi{r') within Eq. 12 



A'^ihi] = 4/ dx'dy'dz'h,{x' -L.j2,y' ~Ly/2,z' -L,/2) 
V Jv 

HxTTx' riy-Ky' n^TTz' 
X cos cos — cos — , (14) 

Lx Ly Lz 

(where the integral is over the range x' G [0, L^], y' £ [0, Ly], z' e [0, L^]), giving 
finally 

CO 

9u{r) = J2 ^'^'^^'ClClCl^ (-4"{/zJ)' , (15) 

^riy .Tiz —0 

with the characteristic time Tc for a given spatial mode defined by 

Putting this expression into Eq. [3] and performing the integral over t gives 

h.{u)= £ ClClCl ^^^i {AHh^})\ (17) 

Substitution into Eqs. [T] and [2] results in complete expressions for the longitu- 
dinal and transverse relaxation times. 

For arbitrary fields including field maps, the cosine-transform amplitudes A 
can be numerically calculated by fast discrete cosine transform (DCT) over hq(f) 
sampled at a sufficient number of points throughout the cell volume. In this 
case the summation over Uq is truncated accordingly, and accuracy is checked 
by increasing the number of sample points and comparing results. 

2.1. Extension to dressed spins with non-uniform dressing field 

An RF magnetic field with amplitude Bi applied transverse to the holding 
field Bq modifies the effective precession frequency of a particle. In terms of 
dimensionless "dressing parameters" [8] 

x^^, y^^, (18) 

in the limit Y 1 the effective gyromagnetic ratio becomes [7] 

j,s = lJo{X). (19) 

Thus for a dressing field with spatially varying amplitude Bi{r) — {Bi)+5Bi{r), 
the equivalent variation SBq in the holding field Bq is given by [T] 

7eff<5Bo ^ Sj,sBo = Bo'^SX = Bo-/Ji{X)6X, (20) 
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The expressions from the previous sections can be used to calculate T2 for the 
dressed spin with a non-uniform dressing field by setting 

7 7off, (21) 

i^Q 7offSo, (22) 

Ki^ ^ r|^<5i?i(f)+Ji?o(rl, (23) 
where the variation 6Bq in the holding field itself has been included in hz- 



3. Linear electric field frequency shift 

A spin moving though an electric field experiences a motional magnetic field 
that may, in conjunction with gradients of the magnetic field, produce a shift 
in the precession frequency dependent on the electric field direction and magni- 
tude. [3] Of particular concern in searches for electric dipole moments are effects 
that are linearly proportional to the electric field E. These may mimic effects 
expected for an electric dipole moment, thereby creating a "false EDM." 

As shown by Lamoreaux and Golub [5] , the linear- in-electric-field frequency 
shift for spins in a confined volume is given by the expression 

1 /•* 

SuJe = / ^'^ COSUJoT{{^x{t)^y{t - t)) - {LO^(t ~ T)ljJy{t))} , (24) 



2 7o 

where the perturbations can be written more generally as 

E E 

^X{t) = lK{t) + 1—Vy{t), UJy{t) ^ jhy{t) - -f—V^{t). (25) 

c c 
Expanding the expression for 6uje and keeping only terms linear in E results in 

2c 

-{hy{t ~ T)Vy{t)) - {Kit ~ T)v,{t))}. (26) 



= -7^ I dr cosujaT{{hy{t)vy{t - t)) + {hx{t)vx{t - t)) 



The cosine-transform method developed in the present work can be used to 
compute Eq. [26] in the diffusion limit for the nEDM cell geometry. While the 
expressions, 

y{t)=yo+ Vy{t')dt', y{t-T)=yo+ vy{t')dt', (27) 

Jo Jo 

were used in Ref. [51 to eliminate y in favor of an expression with Vy, here we 
remove the velocity components from the correlation functions and use instead 
the Fundmental Theorem of Calculus and the above expression for y{t — T). The 
correlation functions can then be written 

{hyit)vyit - r)) = - — {hyit)yit - t)), (28) 
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{hy{t ~ T)Vy{t)) = {hy{t)Vy{t + t)) = ^ {hy{t)y {t + T)) . {29) 

In the latter expression, averages are assumed to be independent of the overall 
time offset i, as appropriate for a stationary problem. 

The expressions in Section [2] are modified to give the correlation function 
in the diffusion limit in terms of the cosine-transform components of hy (r) and 

yif), 

{hy{t)y{t-r))= e-^^^^ClClClAHhyjAHy}. (30) 

{hy{t)vyit-r))= J2 ^e-^'^^ClClClAHhy}AHy} (31) 

Performing the integral in Eq. [26] we have 
f°° 1 

/ dr cos uoT{hy{t)vy{t~T)) = Y T^rTZ^A(^X(^lA{^y\^'{v\^ 

Jo -L + \-^0'r J 



(32) 

leading to an expression for the frequency shift, 



duJE = 



C ^ 1 + {iOoT^Y 

i^y 

(33) 

The summation can be reduced by computing the cosine transform components 
of x{r) and y{r) analytically, 

^ ^ [0 otherwise. ^ ' 

The result for the frequency shift in the diffusion approximation is 



6uj 



„„^3.... 1 + (^ori°'"'"°^)2 „.^3,... "''^^ 1 + (c^ori"-"'"))^ 

(35) 



4. Example application employing the discrete cosine transform 

The bulk of the computational effort required for practical application of the 
present technique is in finding the cosine transform amplitudes A" of the field 
non-uniformities. The form of Eq. [14] is amenable to numerical computation 
with Fast Fourier Transform software libraries. The example below uses the 
multidimensional discrete cosine transform (DCT) feature of the freely-available 
software library FFTW3 10 . Input to the DCT for each field component is an 
array of field perturbations hq{r) sampled over the cell volume at N^NyN^ 
grid points, and the output is an array of DCT amplitudes which, after scaling 
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by l/{8NxNyNz), correspond to the desired amplitudes A^. The summations 
in Eqs. [TT] and [35] are truncated according to the ampHtudes available from the 
DCT. Accuracy of the result can be checked by increasing the number of sample 
points and repeating the computation. 

We calculate the effect of a superconducting rod placed near the cell in an 
otherwise uniform holding field. For an infinite-length rod along the z axis and 
through the origin in a magnetic field Bq applied along the x axis, the net field 
around the rod is pJJ 

Bp = )BoCos0, (36) 

P 

B^ = (1 + — )Bosin(/., (37) 
P^ 

where a is the rod radius and (p, 4>) are polar coordinates in the xy plane. While 
the cosine transforms of B^, By derived from these equations (appropriately 
translated to the desired location of the superconducting rod) could perhaps 
be calculated analytically, here the equations are used to generate a field map 
that is subsequently run through the machinery to produce values for Ti, T2, 
and 5uje- Results are shown in Figure [l] with physical parameters given in the 
caption. 



5. Conclusion 

A method to calculate spin relaxation times and the linear-in-£' frequency 
shift in the diffusion approximation was presented. The technique is based on 
the observation that, for a particle diffusing in a rectangular cell, the correlation 
function of a position-dependent field is the weighted sum over the squared 
cosine-transform components of the field (see Eq. 15 ) ^ 



As the formulation is intended for practical computation in a rectagular cell, 
the result is for the complete 3-dimensional geometry. It could also be applied to 
a cylindrical cell if the field variation in coordinates transverse to the cylinder 
axis can be neglected. Actual computation is done using fast discrete cosine 
transforms of the field components. This method could be used for example 
in magnet coil design optimization: given a field map, Ti, T2, and Suje can be 
quickly evalulated and combined into a figure of merit. 
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•^During preparation of this manuscript, independent work was published based on es- 
sentially the same observation. These authors point out that their result may be used as a 
probe for possible unknown spin interactions and applied it to improve the limits on axion-like 
interactions with the cell walls |13| . 
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Relaxation times vs. distance of rod from cell 
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Frequency shift vs. distance of rod from cell 
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Figure 1: Calculated |(a)| relaxation times Tl, T2 and |(b)| frequency shift Slje due to the 
magnetic field distortion from a 3-mm-diameter, infinitely-long superconducting rod parallel 



to the 2-axis and placed near the cell. The parameters are 7 



-20393.963 radG-^s" 



(gyromagnetic ratio for ^He), D = 428 cm^/s (corresponding to a temperature of 450 mK \12\). 
Bq = 10 mG (applied parallel to the x-axis), E = 100 kV/cm (parallel to Bq), cell dimensions 



: 7.8 cm, Ly 



10 cm, L-, 



■ 40 cm. 
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